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The spontaneous flapping of a flag in a steady flow can be used to power an output circuit using 
piezoelectric elements positioned at its surface. Here, we study numerically the effect of inductive 
circuits on the dynamics of this fluid-solid-electric system and on its energy harvesting efficiency. 

In particular, a destabilization of the system is identified leading to energy harvesting at lower flow 
velocities. Also, a frequency lock-in between the flag and the circuit is shown to significantly enhance 
the system’s harvesting efficiency. These results suggest promising efficiency enhancements of such 
flow energy harvesters through the output circuit optimization. 


I. INTRODUCTION 

Flow induced instabilities and vibrations have recently 
received a renewed attention as potential mechanisms to 
produce electrical energy from geophysical flows (wind, 
tidal currents, river flows, etc.). They indeed enable a 
spontaneous self-sustained motion of a solid body which 
can be used as a generator, effectively converting this 
mechanical energy into electrical form [1-5]. 

A canonical example of such instability is the flapping 
of a flexible plate in an axial flow (e.g. a flag), which 
has been extensively investigated for its rich and com¬ 
plex dynamics [6]. The origin of this instability lies in 
a competition between the destabilizing fluid force and 
the stabilizing structural stiffness. Beyond a critical flow 
velocity U c , the flag becomes unstable, leading to large 
amplitude self-sustained flapping [7 14]. 

Energy harvesting based on flapping plates may follow 
two routes: producing energy either from the displace¬ 
ment [15] or from the deformation of the plate [16-18]. 
The latter has recently been the focus of several studies 
based on active materials [19]. 

Piezoelectric materials, considered in this article, pro¬ 
duce electric charge displacement when strained [20], a 
“direct piezoelectric effect” that effectively qualifies them 
as electric generators. This electric charge can be used in 
an output circuit connected to their electrodes, as in vi¬ 
bration control applications [21]. Piezoelectric materials 
also introduce a feedback coupling of the circuit onto the 
mechanical system: any voltage between the electrodes 
creates an additional structural stress that modifies its 
dynamics (inverse piezoelectric effect). 

The concept of piezoelectric energy generator has re¬ 
ceived an increasing amount of interest in the last 20 
years [22, 23]. Its basic idea is to convert ambient vibra¬ 
tion energy to useful electric energy through piezoelec¬ 
tric materials implemented on vibration sources. Many 
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researchers have contributed to this field in order to im¬ 
prove the efficiency of such energy-harvesting systems 
[24-27]. Some studies show that simple resonant circuits, 
i.e. resistive-inductive circuits combined with the piezo¬ 
electric material’s intrinsic capacitance [20], offer promis¬ 
ing opportunities to achieve high efficiency [28, 29]. 

Flow energy harvesting can be achieved by exploiting 
the unsteady forcing of the vortex wake generated by an 
upstream bluff body to force the deformation of a piezo¬ 
electric membrane [16, 17, 30]. Fluid-solid instabilities 
offer a promising alternative as they are able to generate 
spontaneous and self-sustained structural deformation of 
the piezoelectric structure, e.g. cross-flow instabilities 
[29, 31, 32]. In their work, De Marquis et al. [29] used a 
resistive circuit and a resistive-inductive one, and found 
in addition to the beneficial effect of the resonance to the 
energy harvesting, that a resistive-inductive circuit may 
also affect the stability of the vibration source. How¬ 
ever, the resonant circuit’s influence on the structure’s 
dynamics was not reported in this work. 

Fluid-solid instabilities in axial flows, including the 
aforementioned flapping flag instability, are also stud¬ 
ied in the context of piezoelectric energy-harvesting [33- 
36]. In particular, Michelin & Doare [34, 36] considered a 
piezoelectric flag coupled with a purely resistive output. 
They observed moderate efficiency, which is maximized 
when the characteristic timescale of the circuit is tuned 
to the frequency of the flag. A significant impact of the 
circuit’s properties on the fluid-solid dynamics was also 
identified. 

The present work therefore focuses on the coupling 
of the fluid-solid system, i.e. the flapping piezoelectric 
flag, to a basic resonant circuit (resistive-inductive loop 
with the piezoelectric material’s intrinsic capacitance). 
Resonance is expected when the flapping frequency u> of 
the flag, forcing the circuit, matches the circuit’s natural 
frequency. Using linear stability analysis and nonlinear 
numerical simulations of the fluid-solid-electric coupled 
system, we investigate the impact of such resonance on 
the dynamics and on the amount of energy that can be 
extracted from the device (i.e. the energy dissipated in 
the resistive elements). This explicit description of both 
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the fluid-solid and the electric systems’ dynamics pro¬ 
vides a deeper and more accurate insight into the energy 
harvesting process than its classical modeling as a pure 
damping [4, 15, 18]. 


II. FLUID-SOLID-ELECTRIC MODEL 
A. Fluid-solid coupling 

The coupled system considered here is a cantilevered 
plate of length L and span H , placed in an axial flow 
of density pf and velocity Uf. The flag’s surface is cov¬ 
ered by pairs of piezoelectric patches (Fig. la). Within 
each pair, two patches of reversed polarities are con¬ 
nected through the flag, the remaining electrodes being 
connected to the output circuit [34, 37]. The resulting 
three-layer sandwich plate is of lineic mass p and bend¬ 
ing rigidity B. We restrict here to purely planar de¬ 
formations (bending in the ^-direction and twisting are 
neglected). The flag’s dynamics are described using an 
inextensible Euler-Bernoulli beam model forced by the 
fluid: 


pX = (T t - M! n)' + F f | (1) 

X' = r (2) 

with clamped-free boundary conditions: 

at s = 0:X = X = 0, (3) 

at s = L : T = M = M' = 0. (4) 


Here T and A4 are respectively the tension and the bend¬ 
ing moment. Throughout this article, ' and ' denote 
derivatives with respect to t and s, respectively. The fluid 
loading Ff is computed using a local force model from the 
relative velocity of the flag to the incoming flow: 

U n n + U t t = X — U f e x . (5) 


The present fluid model includes two different contribu¬ 
tions. The first one results from the advection of the 
fluid added momentum by the flow, an inviscid effect, 
and can be obtained analytically in the slender body 
limit through the Large Amplitude Elongated Body The¬ 
ory [38]: 


Freact — TYl apH 


(Jn - ( u n u T y + -uiff 


n. (6) 


Candelier et al [39] recently proposed an analytic proof 
of this result, and successfully compared it to RANS sim¬ 
ulations for fish locomotion problems. These authors also 
stated that in the case of spontaneous flapping, it is nec¬ 
essary to account for the effect of lateral flow separation, 
which is empirically modeled by the following term [40]: 


F 


resist — 


-- 2 pHC d \U n \U n n, 


( 7 ) 



FIG. 1. Flapping piezoelectric flag in a uniform axial flow, 
(a): 3-dimension view; (6): zoom of the circled area in ( x,y ) 
plane; (c): equivalent circuit of a piezoelectric pair connected 
with a parallel RL circuit. 


where Cd is the drag coefficient for a rectangular plate 
in transverse flow. It is important to mention that the 
flapping flag dynamics considered here implies a large 
Reynolds number. A 10 cm long/wide flag in a wind 
flowing at around 5 m/s leads to Re ~ 10 4 , while in wa¬ 
ter we would have Re ~ 10 5 for a flow velocity around 
1 m/s. These values of Re are sufficiently large to jus¬ 
tify a constant value of Cd [41]. A Re-dependence could 
however be introduced to extend the applicability of this 
model to intermediate Re. 

The fluid forcing is the sum of these two terms: 

Ff = Freact T Fresist ■ (8) 

The applicability of this result to flapping flag was con¬ 
firmed experimentally, at least up to an aspect ratio 
H/L = 0.5, the value considered in our work [42]. 


B. Piezoelectric effects 

Piezoelectric patch pairs are positioned on either side 
of the plate, with opposite polarity. This guarantees that, 
during the plate’s bending, the two patches reinforce each 
other, rather than cancel out. Each piezoelectric pair is 
connected to a resistance and an inductance in parallel 
connection (Fig. lc). 

We focus on the limit of continuous coverage by in¬ 
finitesimal piezoelectric pairs [34, 36]. Within this limit, 
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the electric state of the piezoelectric pairs is characterized 
by the local voltage v and lineic charge transfer q , which 
are continuous functions of the streamwise Lagrangian 
coordinate s. The electrical circuits are characterized by 
a lineic conductance g , and a lineic inductive admittance 
l/l. The electrical charge displacement across a piezo¬ 
electric pair resulting from the direct piezoelectric effect 
is given by: 


q = x e' + cv, (9) 

where y is a mechanical/piezoelectric conversion coeffi¬ 
cient and c is the lineic intrinsic capacity of a piezoelectric 
pair [43]. Equation (9) shows that the effect on the circuit 
of the piezoelectric components is that of a current gen¬ 
erator with an internal capacitance (Fig. lc). The charge 
conservation of the resulting RLC circuit leads to: 

v + glv + Iq = 0. (10) 

The inverse piezoelectric effect manifests as an added 
bending moment, so that the total bending moment in 
the structure is given by [37]: 

M = B9' — x v - ( 11 ) 


Finally we define the harvested energy as the time av¬ 
erage of the total rate of dissipation in the resistive ele¬ 
ments, which is formally given as: 


3 ? = 



and the efficiency is defined as: 


V = 



( 12 ) 


(13) 


^■ef is the kinetic energy flux of fluid passing through 
the cross section occupied by the flag: 

^ref = \p U Io XJ^H, (14) 

where s3 is the peak-to-peak amplitude of the flapping 
flag. 

In the following, the problem is nondimensionalized 
using the elastic wave velocity U s = \JB/L 2 p as charac¬ 
teristic velocity. L , L/U s , pHL 2 , U S \Jp/c, U ay f\xc are re¬ 
spectively used as characteristic length, time, mass, volt¬ 
age and lineic charge. As a result, six non-dimensional 
parameters characterize the coupled system: 


p f HL 
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with M* the fluid-solid inertia ratio, U* the reduced flow 
velocity, and H* the aspect ratio. The piezoelectric cou¬ 
pling coefficient, a, characterizes the fraction of the strain 


energy transferred to the circuit, and as such critically 
impacts the energy harvesting performance. Finally (3 
and ui o characterize respectively the resistive and induc¬ 
tive properties of the circuit. 

The effect of the mechanical parameters, M *, U* and 
H *, on the flapping flag dynamics have been extensively 
studied in the literature [34, 36]. In the following, we 
focus specifically on the dynamical properties of the cir¬ 
cuit and maintain M* = 1 and H* = 0.5 throughout 
this study. Unless stated otherwise, we will also con¬ 
sider a = 0.3, a value consistent with existing material 
properties [34]. The effect of varying the piezoelectric 
coupling will also be briefly discussed. The full nonlin¬ 
ear dynamics of the coupled system are now described in 
non-dimensional form by: 

x = M*(Tt)' - (9"n)' + a (v'n)' + M* f { n, (16) 

x' = T, (17) 

j3v + v + (3tOgV + a/39' = 0, (18) 

and the nondimensional boundary conditions are: 

at s = 0:x = x = 0, (19) 

at s = 1 : T = 9' - av = 9” - av' = 0. (20) 


The non-dimensional tension T is computed using the 
inextensibility of the beam [44]. 

Finally the nondimensional fluid loading is obtained 
as: 


ft = ~m a H* 


Un ('U n 'U T ) T yU n 


^ C'd \u n \u n 


( 21 ) 

with m a = 7t/4 and Cd = 1.8, the added mass and drag 
coefficients for a rectangular plate in transverse flow [45, 
46]. 


III. CRITICAL VELOCITY 

The critical velocity U* is defined as the minimum flow 
velocity above which self-sustained flapping can develop 
and energy can be harvested. In this part, the influence 
of the RL loop on U* is studied using linear stability 
analysis, in the limit of small vertical displacement, i.e. 
y <C 1, allowing linearization of Eqs. (16)-(18). The 
resulting linear equations are: 

(1 + M a ) y + 2 M a U*y' + M a U* 2 y" + y"" - av" = 0, 

( 22 ) 

/3v + v + [3 wqV + a(3y" = 0, (23) 

where M a = ttM* H* /4. 

Eqs. (22) and (23) are then projected onto the funda¬ 
mental beam modes and their second derivatives, respec¬ 
tively, and recast as an eigenvalue problem. The coupled 
system is unstable if one of its eigenfrequencies has a 
positive imaginary part. 
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FIG. 2. Evolution of the critical velocity with too at a = 0.3 
and [3 = 0.05 (solid), [3 = 1 (dash-dot), [3 = 4 (dashed), [3 = 8 
(dotted). U* = Uc is plotted (dash-dot, gray) as a reference. 


The evolution of U* with ujq is computed using lin¬ 
ear stability analysis and is shown on Fig. 2. For inter¬ 
mediate values of w 0 (3 < w 0 < 10), we observe a sig¬ 
nificant destabilizing effect of inductance that increases 
with /?, as the circuit becomes dominated by inductive 
effects. For small [3, however, no such destabilization is 
observed, as the inductance plays little role in this re¬ 
sistive limit. These results highlight a major benefit of 
the circuit’s inductive behavior: the instability thresh¬ 
old may be lowered, resulting in energy harvesting with 
slower flow velocity. For u 0 2> 1> the inductance acts 
as a short-circuit, and U* converges, regardless of /3, to 
f7°, the critical flow velocity without coupling (a = 0). 
For w 0 -C 1, the effects of inductance are negligible, and 
U* > [/?, illustrating the stabilizing effect of the resis¬ 
tance. Note that a destabilizing effect of the resistance 
can be observed at higher values of M* [34], and in more 
general cases of damping [47]. 

To determine the origin of this destabilization, Fig. 3 
shows the evolution of the two most unstable pairs of 
eigenvalues with ui o at U* = 10, which is lower than 17° 
but higher than the minimum critical velocity (Fig. 2). 
Starting from uo > 1 and decreasing loq, the electrical 
circuit evolves successively from short circuit, to RLC 
loop, and finally to a purely resistive circuit. Instability 
occurs when the imaginary part of any mode becomes 
positive. 

In the absence of coupling (uj 0 1), both pairs con¬ 

sist of two eigenvalues: (i) an electrical mode with a 
frequency equal to tog; (h) a mechanical mode, with a 
frequency independent of ujg. Decreasing ujq leads to in¬ 
teractions between the electrical and mechanical modes, 
successively within each pair. This interaction destabi¬ 
lizes the mechanical mode, leading to the flag’s instabil¬ 
ity at intermediate uig (Fig. 3). Note that this interaction 
within other pairs also leads to an increase of 3(w) for 
the mechanical mode, but does not lead to instability (at 
least for M* = 1). 
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FIG. 3. (color online) Evolution with cjo of (a) the imaginary 
part (growth rate) and (6) the real part (frequency) of the 
second and third pairs of eigenvalues for a = 0.3, U* = 10 
and [3 = 4. In both (a) and (6), solid lines represent me¬ 
chanical modes and dashed lines represent electrical modes. 
The destabilized pair is plotted in red and the other plotted 
in black. In (6), too is plotted with thin gray solid line as a 
reference. 


IV. NONLINEAR DYNAMICS AND ENERGY 
HARVESTING 

Above the critical velocity, the unstable coupled sys¬ 
tem experiences an exponential growth in its amplitude, 
which eventually saturates due to nonlinear effects. A 
direct integration of the fluid-solid-electric system’s non¬ 
linear equations Eqs. (16)-(18) is performed using an 
implicit second order time-stepping scheme [44]. The 
flag is meshed using Chebysliev-Lobatto nodes, and a 
Chebyshev collocation method is used to compute spatial 
derivatives and integrals. At each time step, the result¬ 
ing nonlinear system is solved using Broyden’s method 
[48]. The simulation is started with a perturbation in 
the flag’s orientation ( 0(s,t = 0) ^ 0), and is carried 
out over a sufficiently long time frame so as to reach a 
permanent regime. 

The reduced flow velocity is chosen at U* = 13, a value 
sufficiently higher than the critical velocity U*. The flag’s 
behavior is observed to drastically differ with varying tog 
(Fig. 4). When tog is within the range of destabiliza¬ 
tion, the flag undergoes a remarkably larger deformation 
(Fig. 46) than with other values of oog (Fig. 4a). 

Figure 5a shows the evolution of the efficiency with uiq, 
and demonstrates that this increased flapping amplitude 
indeed leads to a significant efficiency improvement: a 
maximum efficiency of 6% is obtained here, significantly 
higher than the optimized efficiency obtained at M* = 1 
and U* = 13 without inductance (~ 0.1%) [36]. 

Within the high efficiency range, the flapping fre¬ 
quency is deviated and locks onto the natural frequency 
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FIG. 4. Flapping motion of flags at a = 0.3, j3 = 4, U* = 13 
and (a) u>o = 3.25 (no lock-in), (6): u>o = 4.12 (lock-in). 

of the circuit, wo (Figure 55). A frequency lock-in is 
therefore observed here, similar to the classical lock-in 
observed in Vortex-Induced Vibrations (VIV) [49, 50]: 
the frequency of an active oscillator (the flag) is dictated 
by the natural frequency of a coupled passive oscillator 
(the circuit). The lock-in range is extended by a re¬ 
duction of the circuit’s damping (1//?), consistently with 
what is observed in VIV for varying structural damping 
[51]. The lock-in range leading to high efficiency coin¬ 
cides with the range of wo associated with the desta¬ 
bilization by inductance (Fig. 2). This suggests that a 
coupled piezoelectric flag, once destabilized by inductive 
effects, may flap at a frequency close to the natural fre¬ 
quency of the circuit. As a result, a permanent resonance 
takes place between the flag and the circuit, leading to 
increased flapping amplitude and harvesting efficiency. 

By varying within the lock-in range, Fig. 5 shows 
that when the frequency of the output circuit matches 
the short-circuit natural frequency of the flapping flag 
(w ~ 17.5), the maximal efficiency is obtained for every 
value of f3. This observation highlights again the interest 
of exciting piezoelectric structures at their natural fre¬ 
quencies for energy harvesting, as suggested by previous 
studies, where maximal efficiency is observed when the 
external forcing resonates with the piezoelectric system 
[17, 52, 53]. The existence of a lock-in extends this ef¬ 
fect to a larger range of parameters, by maintaining the 
system at resonance, hence guaranteeing efficient energy 
transfers from the flag to the circuit. 


V. IMPACT OF PIEZOELECTRIC COUPLING 

A decisive factor is the intensity of piezoelectric cou¬ 
pling, characterized by a in this work. It quantifies the 
proportion of the mechanical work transmitted to the cir- 


FIG. 5. (a) Harvesting efficiency r/ and (6) flapping frequency 
u as a function of u>o for a = 0.3, U* = 13, and (3 = 0.05 
(solid), (3 = 1 (dash-dot), f3 = 4 (dashed), /3 = 8 (dotted). 
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FIG. 6. Evolution of the critical velocity with u>o and /3 = 4, 
M* = 1 and a = 0.1 (dash-dot), a = 0.2 (dashed) and a = 0.3 
(solid). 

cuit via piezoelectric effects. In practice, this coupling co¬ 
efficient a, defined in Eq. (15), depends primarily on the 
materials used for the piezoelectric flag. The importance 
of the coupling factor has been reported by many studies 
on energy harvesting by piezoelectric systems [25, 27]. In 
Refs. [34, 36], a dependence of the harvesting efficiency 
to a 2 was identified. 

In the present work, the influence of a is observed both 
in terms of critical velocity (Fig. 6) and the lock-in range 
(Fig. 7). In Fig. 6, lower critical velocities and larger 
destabilization range in terms of are observed with 
increasing a. Fig. 7 shows that the range of frequency 
lock-in also increases with increasing a. The impact of 
the coupling coefficient a on the system’s performance 
is again identified: a strong piezoelectric coupling is de¬ 
sired so that the beneficial effects of a resonant circuit, 
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FIG. 7. Flapping frequency w as a function of u>o for /3 = 4. 
M* = 1, U* = 13 and a = 0.1 (dash-dot), a = 0.2 (dashed) 
and a = 0.3 (solid). 

namely the destabilization and the frequency lock-in, can 
be obtained. 

In practice, a ~ 0.3 could be expected with large scale 
devices [34], and in general a ~ 0.1 is achievable, for 
example using two 10 cm x 10 cm MFC patches glued 
by an epoxy layer with a thickness of 0.1 mm. A potent 
piezoelectric material, leading to a strong coupling, is 
therefore an essential prerequisite to utilize the lock-in 
phenomenon. 

VI. CONCLUSION 

The results presented here provide a critical and 
new insight on the dynamics of a piezoelectric energy¬ 
harvesting flag. First and foremost, they emphasize how 
the fundamental dynamics of the energy harvesting sys¬ 
tem and of the output circuit may strongly impact the 
motion of the structure and its energy harvesting per¬ 
formance. Also, they identify two major performance 
enhancements associated with the resonant behavior of 
the circuit, namely (i) a destabilization of the fluid-solid- 
electric system, leading to spontaneous energy harvesting 
at lower velocity; and (ii) a lock-in of the fluid-solid dy¬ 
namics on the circuit’s fundamental frequency, resulting 
in an extended resonance and a significant increase of the 
harvested energy. 

This lock-in behavior at the heart of both effects above 
is classically observed in VIV where it is also responsible 
for maximum energy harvesting [5]; it is in fact a general 
consequence of the coupling of an unstable fluid-solid sys¬ 
tem to another oscillator’s dynamics. We therefore ex¬ 
pect that the conclusions presented in the present paper 
go beyond the simple inductive-resistive circuit consid¬ 
ered here, and should be applicable to a much larger class 
of resonant systems. Such systems could be other forms 
of electrical output circuits, or other mechanical oscilla¬ 
tors. 



12 14 16 18 20 22 24 

U* 

FIG. 8. (a) Flapping frequency u and ( b ) harvesting efficiency 
ri as a function of U* for a = 0.3, M* — 1, /3 = 3.5 and 
ujo = 0.1 (dotted), ojo = 29 (solid), ojq = 1000 (dashed). 


The critical impact of the coupling coefficient a on the 
system’s performance is also underlined in the present 
work through its strong influence on the destabilization 
range and the lock-in range. The choice of the piezo¬ 
electric materials is therefore essential in the practical 
achievement of high efficiency. 

This lock-in mechanism also plays a critical role in the 
robustness of the energy harvesting process with respect 
to the flow velocity (Fig. 8). Lock-in indeed persists over 
a wide range of U*, effectively acting as a passive con¬ 
trol of the flapping frequency in response to flow veloc¬ 
ity: while ui increases rapidly in the limit of weak fluid- 
solid-electric interactions, lock-in with the output circuit 
maintains u> « ojq and high harvesting efficiency over a 
large range of flow velocity. Such a control of the flap¬ 
ping frequency was shown to be essential for efficiency en¬ 
hancement and robustness (e.g. to delay mode switches 
or frequency changes) [36]. 

These results nonetheless illustrate the fundamental in¬ 
sights and technological opportunities offered by the full 
nonlinear coupling of a passive resonant system (elec¬ 
tric, mechanical or other) to an unstable piezoelectric 
structure for the purpose of energy harvesting. The lock- 
in phenomenon and the enhanced performance demon¬ 
strated by the coupling between the piezoelectric flag and 
a simple resonant circuit open the perspective of applying 
different kinds of resonant systems to energy-harvesting 
piezoelectric system. The choice of an inductive circuit 
is motivated by its simplicity, while in practice, other de¬ 
signs of resonant circuits may present important advan¬ 
tages over the proposed formulation in our work. Mean¬ 
while, complex circuitry, such as propagative [37] and ac¬ 
tive circuits [54], also represent interesting perspectives 
for future work. 
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